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Heavy-ion collisions generate a huge magnetic field of the order of 10 G for the duration of about 0.2 
fm/c. This time may become an order of magnitude longer if the electrical conductivity of quark matter is 
large. We calculate the conductivity in the regime of high density and show that contrary to naive expectations 
it only weakly depends on the magnetic field strength. 



A large body of experimental and theoretical re- 
search conducted for over more than ten years of oper- 
ation of RHIC and the first runs of heavy ion program 
at LHC has led to a revolutionary change in our view 
on the nature and properties of the produced strongly 
interacting matter. A powerful way to investigate the 
nature of a certain substance is to study its response to 
the external perturbations. Few years ago it has been 
comprehended that in heavy ion collisions we have such 
a tool. The produced quark-gluon plasma is subject to a 
super-strong magnetic field (MF) generated by colliding 
ions which have large electric charges and are moving at 
speed close to the speed of light. At the collision mo- 
ment and shortly after it (r < 0.2 fm/c) the MF reaches 
the value \eB\ > rrv\ ~ 10 18 G , i.e. it is of a typical 
QCD scale [T]. In presence of quark matter (QM) the 
life time of the MF may be several times longer provided 
the electrical conductivity (EC) of QM is large enough 
'[2]. The properties of QM including its transport co- 
efficients depend on the location of the system in the 
QCD phase diagram, i.e., on the value of the tempera- 
ture and chemical potential. At zero chemical potential 
and high temperature the EC has been calculated by 
lattice Monte Carlo method yielding significantly differ- 
ent results [3] .Here we consider a reverse regime of high 
density and low or moderate temperature. Qualitatively 
these conditions will be specified below. Physically such 
situation may be realized in neutron stars or in future 
experiments at NICA and FAIR. On the theoretical side 
use can be made of the ideas and methods developed in 
condensed matter physics [J]. 

For a wide class of systems the EC can be decom- 
posed into two contributions. The first one is Boltz- 
mann, or Drude, and it corresponds to semi-classical 



approximation when the mean free path of a parti- 
cle is much larger than any other microscopical length 
scale and the probabilities of a particle interactions are 
added. Quantum effects enter into this contribution 
only via Fermi or Bose distribution functions. A variety 
of quantum effects, such as the interference of trajec- 
tories, quantum phase transitions, quantized vortices, 
fluctuating Cooper pairs, Landau levels in magnetic 
field, etc., define the quantum contribution to the EC. 
In solid state physics this contribution is called "quan- 
tum correction" [3] since it is inversely proportional to 
(kpl), where kp is the Fermi momentum, I is the mean 
free path, and (kpl) ^ 1 in solids. This is not the 
case in dense quark matter where (kpl) may even ap- 
proach unity (we remind that according to Ioffe-Regel 
criterion at kpl = 1 transition to Anderson localization 
phase occurs). The description of the EC in terms of 
the above two contributions is legitimate for highly dis- 
ordered systems [3] in the vicinity of superconducting 
phase transition [5], and in AdS/CFT correspondence 
0. " 

Our starting point is the general expression for the 
EC in terms of Matsubara Green's functions (GF) jU[S]. 
The EC momentum and frequency dependent tensor 
ai m (q,uik), w fe = 2knT is given by: 
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= 0.15. The symbol (...) implies 
the averaging over the disorder. For Drude EC this 
procedure is performed independently for each GF and 
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reduces to the substitution of the standard Matsubara 
frequency e„ = irT(2n + 1) by e n = e n + ^sgn(e n ), 
where r is the momentum relaxation time [5J |3] . For 
the quantum contribution the averaging gives rise to an 
infinite series of the "fan" diagrams [4j [5j [7] yielding the 
result presented at the end of this paper. The GF has 
a standard form [5] 



G(p, e r , 
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where \x is the quark chemical potential. In the regime 
of high density and moderate temperature the transport 
coefficients are dominated by the processes occurring in 
the vicinity of the Fermi surface. Hence the momentum 
integration in pQ) is performed in the following way 



d 3 P _ i <m 
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where £ = (p 2 + to 2 ) 1 / 2 — /i. Calculation of the Drude 
EC reduces to the evaluation of a one-loop diagram de- 
fined by HH1])- First we perform the tr operation over 
Dirac induces, then integration in the complex £-plane, 
and finally the e n summation. Replacing the Matsub- 
ara frequency by the physical ones u>k = — we 
write down the resulting expression for the frequency 
dependent Drude EC 
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where vf = Pf/^- The antiquark contribution is 
damped near the Fermi surface and dropped in Eq. <j5j) . 
If in (J5j> we replace /i — > jj, + m, and take the limit 
|p| <C m, we arrive to the standard Drude formula with 
v = rapp/it 2 and fi being the non-relativistic chemi- 
cal potential. For orientation purposes let us estimate 
<j{lu — > 0) for the following set of parameters: n = 400 
MeV, T = 100 MeV, v P = l,r = 0.8 fm. One gets 
a ~ 0.04 fm -1 . To our knowledge, the EC has not been 
calculated before in this domain of the QCD phase dia- 
gram. Results obtained at /i = and different values of 
T [3] differ from each other by an order of magnitude. 
Our value a ~ 0.04 fm -1 lies within the interval of the 
EC values given in 

An important point is that Eq. ([5]) contains a large 
parameter fi ~ 300 — 500 MeV. As we shall see this leads 
to the stabilization of the EC in MF up to (eB/ '/i)r ~ 1. 
Let the constant MF B be directed along the z-axis. 
The Drude EC can be evaluated either from the one- 
loop diagram with MF entering into the propagators [5] , 



or from the Boltzmann kinetic equation. Diagrammatic 
calculation is more cumbersome, both methods lead to 
the same result which is an anticipated generalization of 
([5]). We present the result leaving the straightforward 
but lengthy derivation for the forthcoming detailed pa- 
per. In line with the symmetry requirements the EC 
along z-axis remains unchanged while the transverse 
one is equal to 



<7_l(w = 0,Q) = 
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where ao = <ju{<jJ — 0) given by (0, and Q — eB/[i. 
For the set of parameters considered above, namely for 
jit = 400 MeV, r = 0.8 fm we have fi 2 r 2 < 1 up to 
eB < 5m 2 . 

Now we give a cursory glance on the quantum part 
a' of the EC. As already mentioned, it includes the in- 
terplay of various quantum phenomena. In the regime 
under consideration the major role is played by the for- 
mation of the fluctuation (or precursor) Cooper pairs 
110) . Again we refer to the forthcoming paper analyz- 
ing diagrams containing propagators of such pairs. The 
order-of- magnitude estimate of quantum EC is \a'\ ~ 
e 2 /ir 3 l, where I — vpT is the mean free path, the sign of 
a' may be either positive, or negative. When a' < and 
large by the absolute value, the system approaches the 
Anderson localization regime. When B ^ the above 
estimate transforms into \a'\ ~ e 2 /tt^Ib, where Ib is the 
magnetic length, Ib — {eB)^ 1 / 2 . Therefore a' is not 
proportional to eB as might be naively expected from 
the fact that the number of modes degenerating in a 
unit transverse area is proportional to B. 

To summarize, we have for the first time evalu- 
ated the EC of dense relativistic quark matter with low 
and moderate temperature. We have demonstrated the 
rigidity of the EC on the magnetic field. The details 
of the calculations will be presented in the forthcoming 
paper. 

B.K. is indebted to A.Varlamov for illuminating dis- 
cussions and to L.Levitov for the remarks. Support from 
RFBR grants 08-02-92496-NTSNIL-a and 10-02093111- 
NTSNIL-a are gratefully acknowledged. 

Remark and references added in v. 3: 

The factor [1 + \{^) 2 ] in © is not a corollary of 
the £-integration.It arises in more general approaches- 
see rrniT^ .The authors thank E.Megias for this remark. 
We also note that our conclusion on the rigidity of the 
EC on the magnetic field was confirmed by a very recent 
calculation |13) in the instanton liquid model. 
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